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Abstract 

We define incidence matrices to be zero-one matrices with no zero 
rows or columns. We are interested in counting incidence matrices 
with a given number of ones, irrespective of the number of rows 
or columns. A classification of incidence matrices is considered for 
which conditions of symmetry by transposition, having no repeated 
rows/columns, or identification by permutation of rows/columns are 
imposed. We find asymptotics and relationships for the number of 
matrices with n ones in some of these classes as n — > oo. 

1 Introduction 

In this paper we address the problem: How many zero-one matrices are there 
with exactly n ones? Note that we do not specify in advance the number 
of rows or columns of the matrices. In order to make the answer finite, we 
assume that no row or column of such a matrix consists entirely of zeros. We 
call such a matrix an incidence matrix. 

Rather than a single problem, there are many different problems here, 
depending on what symmetries and constraints are permitted. In general, 
we define Fij^i{n) to be the number of zero-one matrices with n ones and no 
zero rows or columns, subject to the conditions 

• i = if matrices differing only by a row permutation are identified, and 
i = 1 if not; 



• j = if matrices with two equal rows are forbidden, and j = 1 if not; 

• k = if matrices differing only by a column permutation are identified, 
and k — 1 if not; 

• / = if matrices with two equal columns are forbidden, and I = 1 if 
not. 

The notation is chosen so that Fij^in) is a monotonic increasing function of 
each of the arguments i, j, k, I. 

By transposition, it is clear that F^ujin) = Fijki(n) for all i,j,k,l,n. 
So, of the sixteen different functions defined above, only ten are distinct. 
However, among the problems with k = i and / = j, we may decide that 
matrices which are transposes of each other are identified, leading to four 
further counting problems <&y(n), for i, j G {0, 1}. 

For example, there are four matrices with n = 2, as shown: 




The first has repeated columns and the second has repeated rows. The third 
and fourth are equivalent under row permutations or column permutations, 
while the first and second are equivalent under transposition. 

Table [T] gives some values of these functions. The values of F im (n) 
are taken from the On-Line Encyclopedia of Integer Sequences |Hj, where 
this appears as sequence A101370 and -Foioi(^) appears as sequence A049311, 
while the values of -Foon i n ) an d -Pom {n) are obtained from a formula in Corol- 
lary 3.3 in [9] using MAPLE. Other computations were done with GAP jH]. 
The counting problems can be re-interpreted in various ways: 
Counting hypergraphs by weight Given a hypergraph on the vertex 
set {xi, . . . , x r }, with edges E±, . . . , E s (each a non-empty set of vertices), the 
incidence matrix A = (aij) is the matrix with (i,j) entry 1 if X{ G Ej, and 
otherwise. The weight of the hypergraph is the sum of the cardinalities of 
the edges. Thus -Foioi i n ) is the number of hypergraphs of weight n with no 
isolated vertices, up to isomorphism; and -Fnoi (n) is the number of (vertex)- 
labelled hypergraphs of weight n. Putting k = 1 corresponds to labelling the 
edges, a less usual notion. Moreover, putting / = corresponds to counting 
simple hypergraphs (those without repeated edges). The condition j = is 
less natural in this respect, but corresponds to forbidding "repeated vertices" 
(pairs of vertices which lie in the same edges). 
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Table 1: Some values of the fourteen functions 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


F ooo(n) 


1 


1 


2 


4 


7 


16 








F oio(n) 


1 


1 


3 


11 


40 


174 








F ww {n) 


1 


2 


10 


72 


624 


6522 








F 00 oi(n) 


1 


2 


4 


9 


18 


44 








F on(n) 


1 


2 


7 


28 


134 


729 


4408 


29256 


210710 


F W oi(n) 


1 


2 


6 


20 


73 


315 








F wll (n) 


1 


3 


17 


129 


1227 


14123 








^bioi(w) 


1 


3 


6 


16 


34 


90 


211 


558 


1430 


Fom(n) 


1 


3 


10 


41 


192 


1025 


6087 


39754 


282241 


Fiin(n) 


1 


4 


24 


196 


2016 


24976 


361792 


5997872 


111969552 


$00 (n) 


1 


1 


2 


3 


5 


11 








$io(n) 


1 


2 


8 


44 


340 


3368 








$01 (n) 


1 


2 


4 


10 


20 


50 








$n (n) 


1 


3 


15 


108 


1045 


12639 


181553 


3001997 


55999767 



Counting bipartite graphs by edges Given a zero-one matrix A = 
(Aij), there is a (simple) bipartite graph whose vertices are indexed by the 
rows and columns of A, with an edge from r% to Cj if A^ = 1. The graph 
has a distinguished bipartite block (consisting of the rows). Thus, F 0W i(n) 
and Fim(n) count unlabelled and labelled bipartite graphs with n edges 
and a distinguished bipartite block, respectively (where, in the labelled case, 
we assume that the labels of vertices in the distinguished bipartite block 
come first); $oi(n) counts unlabelled bipartite graphs with n edges and a 
distinguished bipartition. 

Counting pairs of partitions, or binary block designs A block de- 
sign is a set of plots carrying two partitions, the treatment partition and the 
block partition. It is said to be binary if no two distinct points lie in the same 
part of both partitions; that is, if the meet of the two partitions is the par- 
tition into singletons. Thus, F ioi (n) is the number of binary block designs 
with n plots. Putting % — 1 or k — 1 (or both) corresponds to labelling treat- 
ments or blocks (or both). Combinatorialists often forbid "repeated blocks" 
(this corresponds to putting / = 0) although this is not natural from the 
point of view of experimental design. Similarly j — corresponds to for- 



3 



bidding "repeated treatments" . The functions $jj (n) count block designs up 
to duality (interchanging treatments and blocks), without or with treatment 
and block labelling and/or forbidding repeated blocks and treatments. 

Counting orbits of certain permutation groups A permutation 
group G on a set X is oligomorphic if the number F*(G) of orbits of G 
on X n is finite for all n. Equivalently, the number F n (G) of orbits on or- 
dered n-tuples of distinct elements is finite, and the number f n (G) of orbits 
on n-element subsets of X is finite, for all n. These numbers satisfy various 
conditions, including the following: 

n 

• F*{G) = S(n, k)F k (G), where S(n,k) are Stirling numbers of the 

k=l 

second kind; 

• f n {G) < F n (G) < n\f n (G), where the right-hand bound is attained if 
and only if the group induced on a finite set by its setwise stabiliser is 
trivial. 

For example, let S be the symmetric group on an infinite set X, and A 
the group of all order-preserving permutations of the rational numbers. Then 
F n (S) = f n (S) = f n (A) = 1 and F n (A) = n\ . 

Now if H and K are permutation groups on sets X and Y, then the direct 
product H x K acts coordinatewise on the Cartesian product 1x7. It is 
easy to see that F*(H x K) = F*{H)F*{K). 

Let (xi,yi), . . . , (x n ,y n ) be n distinct elements of XxY. If both X and Y 
are ordered, then the set of n pairs can be described by a matrix with n ones 
in these positions, where the rows and columns of the matrix are indexed by 
the sets {x±, . . . , x n } and {yi, . . . , y n } respectively (in the appropriate order). 
Moreover, if X is not ordered, then we can represent the set of pairs as the 
equivalence class of this matrix under row permutations, and similarly for 
columns. Thus 

F 0W i(n) = f„(S x S), F n oi(n) = f n (S x A), F 1111 (n) = f n (A x A). 

Moreover, the wreath product H wr C 2 is the permutation group on X 2 
generated by H x H together with the permutation r : (xi,x 2 ) i— > (x 2 ,Xi). 
The effect of r is to transpose the matrix representing an orbit. So 

$01 W = /n(SwrC 2 ), $ n (n) = / n (Awr C 2 ). 



4 



Discussion of this "product action" can be found in [Sj and [T2] . 
It is not clear how forbidding repeated rows or columns can be included 
in this interpretation. 



2 The asymptotics of F im (n 

We will use both F(n) and Fim(n) to denote the number of incidence 
matrices with n ones. This is the largest of our fourteen functions, so its 
value gives an upper bound for all the others. Indeed, we will see later that 
Fijki{n) = o(F nn (n)) for k,l) ^ (1, 1, 1, 1). 

It is possible to compute this function explicitly. For fixed n, let my be 
the number of i x j matrices with n ones (and no zero rows or columns). We 
set m o(0) = 1 and F(0) = 1. Then 

55©CM"> 

so by Mobius inversion, 

'k\ (l\ f ij 

i-u m 

i<k j<l 

and then 

F U n(n) = ^2^2 mij . (3) 

i<n j<n 

For sequence a n , b n , we use the notation a n ~ b n to mean lim^oo a n /b n = 
1. It is clear from the argument above that 

Fmi(n) < ( n2 ) ~ -ft=(pe) n , 



and of course considering permutation matrices shows that 

Fnn(n) > n\ ~ V2im {-) . 



Theorem 2.1 



Fnnin) ~ V' 00850 ' 



4 (log 2) 



2n.+2 
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We remark that for n = 10, the asymptotic expression is about 2.5% less 
than the actual value of 2324081728. 

We have three different proofs of Theorem 12.11 One proof will be given 
in its entirety and the other two will be briefly sketched. Their full details 
can be found in 0. We use the method of the first proof to bound Fn i in 
Section 5. The ideas behind the third proof lead to a random algorithm for 
generating incidence matrices counted by Fim(n) and by $n(n). The ran- 
dom algorithm provides an independent proof of the expression for Fim(n) 
used in the first proof. 



First proof This proof uses a procedure which, when successful, gener- 
ates an incidence matrix uniformly at random from all incidence matrices. 
The probability of success can be estimated and the asymptotic formula for 
Fiui(ri) results. 

Let R be a binary relation on a set X. We say R is reflexive if (x, x) G R 
for all x G X. We say R is transitive if (x, y) G R and (y, z) G R implies 
(x, z) G R. A partial preorder is a relation R on X which is reflexive and 
transitive. A relation R is said to satisfy trichotomy if, for any x, y G X, 
one of the cases (x, y) G R, x = y, or (y, x) G R holds. We say that R is a 
preorder if it is a partial preorder that satisfies trichotomy. The members of 
X are said to be the elements of the preorder. 

A relation R is antisymmetric if, whenever (x, y) G R and (y, x) G R 
both hold, then x = y. A relation R on X is a partial order if it is reflexive, 
transitive, and antisymmetric. A relation is a total order, if it is a partial 
order which satisfies trichotomy. Given a partial preorder R on X, define a 
new relation S on X by the rule that (x, y) G S if and only if both (x, y) and 
(y, x) belong to R. Then S is an equivalence relation. Moreover, R induces 
a partial order R on the set of equivalence classes of S in a natural way: if 
(x, y) G R, then (x, y) G R, where x is the S'-equivalence class containing 
x and similarly for y. We will call an S'-equivalence class a block. If R is a 
preorder, then the relation R on the equivalence classes of S is a total order. 
See Section 3.8 and question 19 of Section 3.13 in jlj for more on the above 
definitions and results. Random preorders are considered in jjj. 

Given a preorder on elements [n] := {1,2, ...,n} with K blocks, let 
Bi, B 2 , . . . , Bfc denote the blocks of the preorder. Generate two random 
preorders uniformly at random B\, B%, . . . , Bk and B[, B' 2 , . . . , B' L . For each 
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1 < 2 < j < n, define the event Dij to be 

Di j = {for each of the two preorders % and j are in the same block}. 
Furthermore, define 

l<i<j<n 

where the indicator random variables are defined by 

j f 1 if D it j occurs; 

1,3 1 otherwise. 

If W = 0, then the procedure is successful, in which case B k R -B[ consists of 
either or 1 elements for each 1 < k < K and 1 < I < L. If the procedure 
is successful, then we define the corresponding K x L incidence matrix A by 



A 



k.i 



1 ifB fc n£;^0; 

o if B k n b; = 0. 



It is easy to check that the above definition of A in fact produces an incidence 
matrix and that each incidence matrix occurs in n\ different ways by the 
construction. It follows that 

_ PjnfnW = 0) 

fini{n) j , 

nl 

where P(n) is the number of preorders on n elements if n > 1 and -P(O) = 1. 

It is known (see jT], for example) that the exponential generating function 
of P(n) is 

^ n! 2-e 2 w 

n=0 

The preceding equality implies that P(n) has asymptotics given by 

n+l 

2 V.log2 



(5) 



where, given sequences a n , b n the notation a n ~ b n means that lim^^oo a n /b Ti 
1. It remains to find the asymptotics of F(W = 0). 
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The rth falling moment of W is 
E(W) r = EW(W -!)■■■(]¥ -r + l) 



= E E W"ArJr (6) 



= E £ 7^ ■ • ■ 7 VJV + E QT* 7 njl ■ • ■ 4., ir )(,7) 




.all i s and j' s different 



pairs (isjs) different 




with \ defined to be the sum with all pairs (i s ,j s ) different, but not all 



i s , j s different. 

First we find the asymptotics of the first term in (JJJ). For given sequences 
ii, i 2 , . . . , i r , 31,32, ■ ■ ■ ,jr, the expectation E(/j 1J1 • • • 7 irJr ) is the number of 
ways of forming two preorders on the set of elements [n] \ {ji,j2, ■ ■ ■ ,j r } 
and then for each s adding the element j s to the block containing i s in both 
preorders (which ensures that D is j a occurs for each s) and dividing the result 
by P{n) 2 . Since the number of ways of choosing i l5 i 2 , . . . ,i r , ji,j2, ■ ■ ■ ,j r 
equals A , This gives 



where we have used ©• 

The second term is bounded in the following way For each sequence 
(hi (*2j 32), ■ ■ ■ 1 {isijs) i n the second term we form the graph G on vertices 
Us=i{*«> 3s} with edges {J r s=1 {{i s , js}}- Consider the unlabelled graph G' 
corresponding to G consisting of v vertices and c components. The number 
of ways of labelling G' to form G is bounded by n v . The number of preorders 
corresponding to this labelling is P(n — v + c) because we form a preorder on 
n — v + c vertices after which the vertices in the connected component of G 
containing a particular vertex get added to that block. Therefore, we have 





G' 



£0 („»"), 



G' 
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where the constant in O (n 2c ~ v ) is uniform over all G' because v < 2r. Since 
at least one vertex is adjacent to more than one edge, the graph G is not a 
perfect matching. Furthermore, each component of G contains at least two 
vertices. It follows that 2c < v and, as a result, 

E(Y,*I ilA ...I ir j r )=0(n- 1 ). 
The preceding analysis shows that 

(log2) 2 ^ 



E(W) r < 

for each r > 0. The method of moments implies that the distribution con- 
verges weakly to the distribution of a Poisson((log 2) 2 /2) distributed random 
variable and therefore 

W>(W = 0) ~ exp (-Q^-^j . (8) 

□ 



Second proof (Sketch) First, the following expression for F im (n) is given 
in terms of the number of preorders on k elements as an alternating sum 
different from and simpler than (J2J): 

1 n 

F ull (n) = -Y / s(n,k)P(k) 2 , 
n! L — ' 

k=l 

where and s(n, k) and S(n, k) are Stirling numbers of the first and second 
kind respectively. As in the first proof, the number of pairs of preorders for 
which the meets of the blocks form a given fc-partition of [n] is k\F(k), so 

n 

P{nf = ^S{n,k)k\F{k), 

k=l 

and we obtain the result by inversion. Next, P(k) is replaced by its asymp- 
totic expression (J3J) with negligible error. Let 

F\n) = \.-J2<n,k)(k\) 2 c k+1 , 

k=l 
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where c = l/(log2) 2 is as in the statement of the theorem. As we have 
argued, F(n) ~ F'{n). 

Now, (— l) n ~ k s(n, k) is the number of permutations in the symmetric 
group S n which have k cycles. So we can write the formula for F'(n) as 
a sum over S n , where the term corresponding to a permutation with k cy- 
cles is (— l) n ~ h (k\) 2 c k+1 . In particular, the identity permutation gives us a 
contribution 



To show that F'{n) ~ Cg{n) as n — > oo, where C = exp(— (log2) 2 /2), we 
write F'(n) = F[(n) + F^n) + F^(n), where the three terms are sums over 
the following permutations: 

F[: all involutions (permutations with a 2 = 1); 

F' 2 : the remaining permutations with k > \n/2~]\ 

F^\ the rest of S n . 

A further argument shows that F[(n) ~ Cg(n), while F^n) = o(g(n)) and 
Fi{n)=o{g{n)). □ 

Third proof (Sketch) If one is interested in asymptotic enumeration of 
F(n), the formula (J2J), being a double sum over terms of alternating sign, 
is on first sight rather unsuitable for an asymptotic analysis. We present 
a derivation of the asymptotic form of F(n) based on the following elegant 
and elementary identity. (This identity and equation J2J) were also derived 

Proposition 2.2 






Proof Insert 




(10) 



into @ and resum using (JTJ). 



□ 
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The sum in (J§J) is dominated by terms where kl 3> n. In this regime, 
using 



-g 2kl 



n J n\ 

and approximating the sum in (jH)) by an integral (cf. Euler-Maclaurin) leads 
to 



F(n) ~ ^Jdkjdl 



(MY 

-g 2fcl 



d/t / dA e n< - logK_Klog2 ^e n( - log A_Alog2 ^e~2^ 



n 2n+2 

4n! 

For n large, the integrals are dominated by a small neighborhood around their 
respective saddles. As e~2^x is independent of n, we can treat the integrals 
separately. Using w(k) = logn — ftlog2, the saddle n s = * s determined 
from w'(k s ) = (X s = analogously). Approximating the integrals by a 
Gaussian around the saddle point gives 



n 2n+2 „„,. ^ / 2tt ^ / 2tt l_ 

2k s X 



F(n) ~ - e nw ( K °) . / e nw(X s ) I — e - 

V) 4n! V«K(«-)r V™K( A s)l 



W 2?1+2 / 3n (i g log 2_l) / 2?r \ p -i(log2) 2 



2 

e 2* 



4n! \ U ra(log2) 2 

which simplifies to the desired result. □ 

3 Generating random incidence matrices 

It is easily shown that (j3J) implies that 

°° k n 

k=0 

Hence, the distribution irk on the natural numbers defined by 

k" 



P(n)2 k + l 



is a probability distribution. The following way of generating preorders uni- 
formly at random was given in [IT]. 
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Theorem 3.1 (Maassen, Bezembinder) Let A be a set of n elements, 
n > 1 . Let a random preorder R be generated by the following algorithm: 

(i) Draw an integer-valued random variable K according to the probability 
distribution ir k . 

(ii) To each a G A assign a random score X a according to the uniform 
distribution on {1,2,..., K}. 

(Hi) Put aRb if and only if X a < X b . 

Then all of the P{n) possible preorders on A are obtained with the same 
probability P{n). 

Incidence matrices counted by -F im (n) can be generated uniformly at 
random by a similar algorithm. Define a integer valued joint probability 
distribution function p^j by 



Theorem 3.2 The following algorithm generates a random incidence matrix 
counted by Fuu(n). 

(i) Draw integer-valued random variables K and L according to the joint 
probability distribution p k \. 

(ii) Choose a 0-1 matrix with K rows, L columns, n l's and KL — n O's 
uniformly at random. 

(Hi) Delete all rows and columns for which all entries are 0. 

Proof Denote a 0-1 matrix with k rows, I columns, and n l's a (k, /)-matrix. 
Denote an incidence matrix with i rows, j columns, and n l's a (i, j)-incidence 
matrix. Now, every (k, Z)-matrix is generated with equal probability 




Pk,l 



Fiin(n) 
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and every (i, j)-incidence matrix is generated from (J (.) [k, /)-matrices. Av- 
eraging over the probability distribution, it follows that every (i, j)-incidence 
matrix is generated with probability 



p(i,j) = 




Using (jTUj) . this sum simplifies to p(i,j) = 1/F im (n). □ 



4 Counting symmetric matrices 

In this section we find the asymptotics for $n(n) and show: 
Proposition 4.1 $n(n) ~ 

Proof Clearly we have 3>n(n) = |(-Fim(n) + S'n(n)), where Sn(ra) is the 
number of symmetric matrices with n ones having no zero rows or columns, 
where repeated rows or columns are allowed and row or column permutations 
are not permitted. So it suffices to show that Su{n) = o(Fmi(n)). 

Now let I(n) be the number of solutions of a 2 = 1 in the symmetric group 
S n . Then we have 

I{n) < S n {n) < I{n)P(n)/n\. 

The lower bound is clear by considering symmetric permutation matrices. 
For the upper bound, our analysis of F im (n) shows that n\Su{n) is the 
number of pairs (Ri, R2) of preorders on {1, . . . , n} such that no two points i 
and j lie in the same block for both preorders, and additionally such that R\ 
and R2 are interchanged by some involution a of {1, ... ,71} (corresponding 
to transposition of the matrix). So instead of choosing R 1 and R 2 , we can 
choose Ri and a and let R2 = R°; there are P(n)I(n) choices, and this is an 
overcount because of the extra condition that must hold on {Ri^R^). 
Now J(n) is just a little larger than \fn\: in fact, 

j( n ) „ —-^ 

A /2e n / 2 "v / ^+ 1 /4 

(see |3 p. 347]). We have seen that P(n)/n\ ~ A(l/log2) n . So the conclusion 
follows from Theorem 12.11 □ 
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It is possible to show that the upper bound for Su(n) is correct, apart 
from a constant factor: 

Proposition 4.2 S n (n) ~ C s ■ I{n)P{n)/n\, where C s = ± e -( log2 ) 2 / 4 « 
0.44341. In other words, if we choose randomly a preorder R and an involu- 
tion a G S n , the probability that no two points lie in the same part in both R 
and R a tends to C s as n — > oo. 

Proof Let /ij = /Jj(rt) be the number of i x % symmetric incidence matrices 
with n ones. Let Sk be the number of k x k symmetric matrices with n ones, 
given by 

•■=g(?)(.-',) : 

where j represents the number of ones off of the diagonal. Then 

k 



and 



i=l 



n oo 

8=1 fc = l 



by (fTUj) . leading to 

^(«) = E^tE(?)( 

fc=0 J=0 v J 7 v 

To compute this sum asymptotically, we approximate for m ^> I ^> 1 



k 

n - 2j 



I J I! " 



The sums are dominated near k ~ nj log 2 and j ~ (n — v / n)/2, so that we 
can justify replacing the binomial coefficients by this approximation. We get 



oo 1 oo (k\3 , n —2j 

_ V A2i_— _ e -i 2 /2©-(n-2j)V2fc 



~ EpE ji (n _ 2j -)! 
fc=o i=o v Jy 



k 



n 



V V 1 r -.?7fc-.7 2 /fc 2 -(n-2.7) 2 /2fc 



fc=0 j=0 J v Jy 
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where in the last step we also replaced (1 — l/k)^ ~ e~^ k . Due to the 
concentration of the sum near k ~ nj log 2 and j ~ (n — y / n)/2, the argument 
of the exponential can be replaced by 

e -j/k-f/k 2 -(n-2jf/2k ^ e - log 2- (log 2/2) 2 _ l e -(log2) 2 /4 _ 

where C s is as in the Proposition. Identifying 

oo 

P{n) = £ 



v k n 



and 



we arrive at 



2 k+1 

k=0 



o ./!(// -2./) !2-' 



S n {n)~C s -I{n)P{n)/n\ 



□ 



One may also generate matrices from S\\ uniformly at random. 
Define an integer valued probability distribution function x/jf. by 

s fc 2-*- 1 



Theorem 4.3 The following algorithm generates a random incidence matrix 
counted by Sii(n). 

(%) Draw integer-valued random variables K according to the probability 
distribution ipk ■ 

(ii) Choose a K x K symmetric zero- one matrix with n ones and K 2 — n 
zeros uniformly at random. 

(Hi) Delete all rows and columns for which all entries are zero. 

The proof of Theorem 14. HI is similar to the proof of Theorem IH.2I 

In general, we have $y(n) = \{Fijij{n) + SV,(n)), where 
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• if % = 1, then Sij(n) is the number of symmetric matrices with n ones 
and no zero rows, where repeated rows are forbidden if j = and 
permitted if j — 1; 

• if % — 0, then Sij(n) is the number of classes of matrices with n ones and 
no zero rows (up to row and column permutations) which are closed 
under transposition, with the same interpretation of j as in the other 
case. 

We do not yet have asymptotics for these. It seems likely that, in all four 
cases, Sij(n) = o(Fijij(n)), so that n) ~ ^Fijij{n). Table |21 gives some 
values of these functions. 



Table 2: Some counts for symmetric matrices and classes 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Soo(n) 


1 


1 


2 


2 


3 


6 










S 01 (n) 


1 


1 


2 


4 


6 


10 










Sio(n) 


1 


2 


6 


16 


56 


214 


866 


3796 


17468 




Su(n) 


1 


2 


6 


20 


74 


302 


1314 


6122 


29982 


154718 



5 The function -FioiifVi) 

Recall that the number of incidence matrices with n ones, no repeated rows 
and matrices equal by row or column permutations unidentified is denoted 
by Fion(n). In this section we will show 

Theorem 5.1 We have 

F 10 n(n) = o(F ini (n)) . 

Proof We will use the probabilistic method and the notation used in the 
proof of Theorem l2.ll The idea behind the proof is to show that the probabil- 
ity tends to that a randomly chosen incidence matrix counted by i 71 im (ra) 
does not have two rows with each containing all zeroes except for a single 
one in the same column. 
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Define Eij, 1 < i,j < n, to be the event that both {i} and {j} are 
blocks in the first preorder and that i and j belong to the same block of the 
second preorder. When W — 0, corresponds to the event that the rows 
corresponding to the blocks containing i and j in the incidence matrix are 
different and contain unique ones appearing in the same column. 

Let P(n, k) be the number of preorders on n elements with k blocks. 
Given a power series f(z) = Y^n=o fn zTl i define [z n ]f(z) = f n . We find that 
for any 1 < i < j < n, 



nEij) = g^l_M ( , + 2)( , + 1) .^zi) 

^ P(n) P{n) 

P( T \ 

= -TT-^^ k + 2 ^ k + l ) p ^- 2 ^ ( 12 ) 

^ > k=l 



Using Lemma 1.1 of jjj, we find that 



71 — 2 / oo 

J2(k + 2){k + l)P(n — 2, k) = {n- 2)\[z n - 2 } J^(n + 2){n + l){e z - 1) 

fe=l \n=0 

d 2 ( U 2 



(n-2)\[z 



n-2i 



du 2 V 1 — u 



u=e z — 1 



When singularity analysis (see Section 11 of [13J) can be applied, as in this 
case, the asymptotics of the coefficients of a generating function are deter- 
mined by the degree of its pole of smallest modulus. Therefore, 

+ 2){k + l)P(n - 2, k) ~ (n - 2)![^ 2 ] ^ 1 - j . 

The singularity of smallest modulus of (2 — e 2 ) -1 occurs at z = log 2 with 
residue 

/^-iog2\ .. /n i 

lim = hm 



2^iog2 y 2 — e z / z^iog2\— e z J 2 
by l'Hopital's rule. Hence, 

n-2 



J2(k + 2)(k+ l)P(n - 2, *) ~ ^-^V~ 2 ](log2 - 



fc=l 
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from which singularity analysis and (jSJ) give 

n-2 



+ 2)(k + l)P(n - 2, k) ~ (n 2)1 (log 2)^ n - 1 n 2 . 



(13) 



fc=i 



The result of using (|T3*j) in (|T2j) is 



log 2 

An 



(14) 



Define X to be 

l<i<j'<n 

so that, conditional on the event {P^ = 0}, the event {X > 0} implies that 
the incidence matrix produced by the algorithm has repeated rows. The 
expectation of X is 



E(X)=(^P(^)~!^n 



We will next show that 

E((W) r nEij) 

The analog of (jOJ is 



(log2) 2 Vlog2 



4n 



(15) 



(16) 



/ 



1,3 1 



E 



E 



pairs (i s ,js ) different 



all i s and j s different 



\ 

i 

J 
\ 

i 

J 



with defined to be the sum with all pairs (i s ,j s ) different, but not all 

i s ,j s different. 
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The first term corresponds to two preorders formed in the following way. 
The i s , j s are first selected. One preorder is formed from the set of elements 
[n] \ {i,j,ji,j2, ■ ■ ■ ,jr}, the element j s is added to the block containing i s for 
each s, and then blocks {i} and {j} are inserted in the preorder. Another 
preorder is formed from the set of elements [n] \ {j, ji, j'2, . . . , j r }, the element 
j s is added to the block containing i s for each s, and then the element j is 
added to the block containing i. As a result, 



E 



/ \ 



pairs (i s ,j s ) different 



"v' P(n~2- r, k) P(n-l-r) 
2»-(n-2-2r)! P(n) l* + W + A J p{n) 



(log 2)^ V log 2 
2 ) An 

where we have used © and the asymptotic form (JT3J). 

The second term is bounded using the same method that was used to 
bound in the proof of Theorem 12. 11 Letting G" be an index over graphs 
on n vertices with two labelled disconnected vertices i and j and n — 2 
unlabelled vertices which is not a matching on the unlabelled vertices, we 
have 



G" fc 



G" 

Ofn" 2 ). 



Consequently we have shown (jlfijl . 

The asymptotics (|14|) and and method of moments argument giving 
(jBJ) imply that 



P(H/ = 0|E w )~exp(-(l2|5!) 
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and therefore an application of Bayes' Theorem with (JHJ) and (|14|) results in 

log 2 

F(Ei j \W = 0) ~ -p-. 

4n 

The observations above result in 

E(X \W = 0)~ ^|^n. (17) 
8 

Comparison of f)15|) and (jl7|) makes it clear that that conditioning on the 
event {W = 0} does not asymptotically affect the expectation of X. 

In a similar way we can find the asymptotics of the conditional second 
falling moment K(X(X — 1) | Y — 0). The unconditioned second moment 
equals 

E(X(X-1)) = F(E h!ji nE i2!j2 ) 

nE iujl nE i2 , h )+ J2 nEn, n nE t2 , 2 ) 

{uji}n{i 2 ,i2}=0 |{nJi}n{j 2 ,i2}|=i 

n\ ^ P(n-A,k) „ Wl Wl Wl ,P(n-2) 



/ 3 gpr^) 



(A; + 1)(A; + 2)(A; + 3) P( " " ! 



P(n) 

An application of singularity analysis as used to derive (|T3|) produces 

E(I(I-l)) = ^ + 0(n). 
Arguing as we did for K(X | W = 0) shows that 

E(A(X - 1) | W = 0) = 2^l!n 2 + O(n); 

64 

we omit the details. 

The variance of X conditioned on W = is 

Var(A \ W = 0) = E(X(X - 1) | W = 0) + E(A | W = 0) - (E(X | W = 0)) 
= o(n 2 ). (18 



2 
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Chebyshev's inequality applied with (f*T7j) and (JTSJ) now gives 

¥(X = | W = 0) = o(l). 

Hence, an asymptotically insignificant fraction of incidence matrices do not 
have repeated rows which implies that iq n(n) = o (-Fini(n)). □ 

6 The functions Faou(n) and Fom(n) 

The function -Fom(n) counts vertex-labelled hypergraphs on n vertices, while 
-^oon (n) counts the simple vertex-labelled hypergraphs. For completeness, we 
include the formulae from the work of Martin Klazar [§| . 

Theorem 6.1 (a) For all n, we have 

w») - EEri(2)E(-^(7). 

Ahn j=i i=l v 1 / m=j v J ' 

Ahn j=l i=l "•— .1 

where A = \ ai 2 a2 ■ ■ ■ l ai is a partition of n with ai > 0. 
(b) For all n, we have 

Foou (n) < Fom (n) < 2F 0011 (n) . □ 

Part (b) raises the question of whether i ? ooii("')/Foiii(n) tends to a limit 
as n — > oo, and particular, whether the limit is 1 (that is, whether almost all 
labelled hypergraphs are simple). 

The paper [5j also gives recurrence relations for the two functions. Klazar 
subsequently showed [TU] that both functions are asymptotically 

(l/log2 + o(l)) n 6(n), 

where b(n) is the nth Bell number (the number of partitions of {1, . . . ,n}). 
Details of the asymptotics of b(n) can be found in [TH]. In particular, since 
b{n)/n\ = (1 + o(l)) n , we see that F m(n) = o(F im (n)), and in fact 

F om (n) = (\og2 + o(l)) n F nil (n). 
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This and the result of the last section, together with the facts that F 1101 (n) = 
Fom(n) and Fmo(n) = F W u(n) and that F^in) is monotone increasing 
in each of i,j,k,l, justify our earlier claim that Fijki(n) = o(Fim(rt)) for 
(i,j,k,l)^ (1,1,1,1). 

7 A rough lower bound for Foooi(^) 

The number Foioi of unlabelled hypergraphs with weight n is not smaller 
than the number of graphs with n/2 edges and no isolated vertices. We show 
that this number grows faster than exponentially. In fact, our argument 
applies to F 0001 , since we use simple graphs. 

Consider simple graphs with m vertices and n edges, where m = o(n) and 
n = o(m 2 ). The number of such graphs, up to isomorphism, is at least 

( m(m ; 1)/2 ) (cm 2 )" c n m 2n 
ml n\ m\ n n m m 

for some constant c. Put F equal to the logarithm of the right-hand side: 

F = c'n + 2n log m — n log n — m log m 
for some constant d . Putting m = c"n/ logn, for some constant c", we get 

F = n log n — 2n log log n + 0(n). 

We conclude: 
Proposition 7.1 For any e > 0, we have 

FoooiH - (a^)^) 

for n > no(e). 

Remark The asymptotics of the number of graphs with no isolated ver- 
tices, having a given number of vertices and edges, has a long history: see P2>j 
for an early paper on this topic, and [2] for a recent result. 
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